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The principle of the Karlsruhe dynamo experiment is closely related to that of the Roberts dynamo 
working with a simple fluid flow which is, with respect to proper Cartesian co-ordinates x, y and z, 
periodic in x and y and independent of 2. A modified Roberts dynamo problem is considered with 
a flow more similar to that in the experimental device. Solutions are calculated numerically, and on 
this basis an estimate of the excitation condition of the experimental dynamo is given. The modified 
Roberts dynamo problem is also considered in the framework of the mean-field dynamo theory, in 
which the crucial induction effect of the fluid motion is an anisotropic a-effect. Numerical results 
are given for the dependence of the mean-field coefficients on the fluid flow rates. The excitation 
condition of the dynamo is also discussed within this framework. The behavior of the dynamo in 
the nonlinear regime, i.e. with backreaction of the magnetic field on the fluid flow, depends on the 
effect of the Lorentz force on the flow rates. The quantities determining this effect are calculated 
numerically. The results for the mean-field coefficients and the quantities describing the backreaction 
provide corrections to earlier results, which were obtained under simplifying assumptions. 
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I. INTRODUCTION 

In the Forschungszentrum Karlsruhe U. Miiller and 
R. Stieglitz have set up an experimental device for the 
demonstration and investigation of a homogeneous dy- 
namo as it is expected in the Earth's interior or in cosmic 
bodies pi. The experiment ran first time successfully in 
December 1999, and since then several series of measure- 
ments have been carried out M 0, 0, |j. It is the second 
realization of a homogeneous dynamo in the laboratory. 
Its first run followed only a few weeks after that of the 
Riga dynamo experiment, working with a somewhat dif- 
ferent principle, which was pushed forward by A. Gailitis, 
O. Lielausis and co-workers (|, [7|. 

The basic idea of the Karlsruhe experiment was pro- 
posed in 1975 by F. H. Busse @, Hj. It is very simi- 
lar to an idea discussed already in 1967 by A. Gailitis 
|jTof . The essential piece of the experimental device, the 
dynamo module, is a cylindrical container as shown in 
Fig. |], with both radius and height somewhat less than 
lm, through which liquid sodium is driven by external 
pumps. By means of a system of channels with con- 
ducting walls, constituting 52 "spin generators", helical 
motions are organized. The flow pattern resembles one of 
those considered in the theoretical work of G. O. Roberts 
in 1972 [jll). This kind of Roberts flow, which proved to 
be capable of dynamo action, is sketched in Fig. ^. In 
a proper Cartesian co-ordinate system [x, y, z) it is pe- 
riodic in x and y with the same period length, which 
we call here 2a, but independent of z. The x and y— 
components of the velocity can be described by a stream 
function proportional to sin(7ra;/a) sin(7ry/a), and the z— 
component is simply proportional to sin(7ra;/a) sin(7n//a). 



When speaking of a "cell" of the flow we mean a unit like 
that given by < x, y < a. Clearly the velocity is contin- 
uous everywhere, and at least the x and y-components 
do not vanish at the margins of the cells. The real flow in 
the spin generators deviates from the Roberts flow in the 
way indicated in Fig. |j| In each cell there are a central 
channel and a helical channel around it. In the simplest 
approximation the fluid moves rigidly in each of these 
channels, and it is at rest outside the channels. We re- 
late the word "spin generator flow" in the following to 
this simple flow. In contrast to the Roberts flow the spin 
generator flow shows discontinuities and vanishes at the 
margins of the cells. 

The theory of the dynamo effect in the Karlsruhe de- 
vice has been widely elaborated. Both direct numer- 
ical solutions of the induction equation for the mag- 
netic field g|, [D| ^ [D| fl6L 0, |l|] as well as mean- 
field theory and solutions of the corresponding equa- 
tions Jl|, |0[ [H], ||, H H| have been employed. We 
focus our attention here on this mean-field approach. 
In this context mean fields are understood as averages 
over areas in planes perpendicular to the axis of the dy- 
namo module covering the cross-sections of several cells. 
The crucial induction effect of the fluid motion is then, 
with respect to the mean magnetic field, described as 
an anisotropic a effect. The a coefficient and related 
quantities have first been calculated for the Roberts flow 
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|22j, |23|, |25j. In the calculations with the spin 
generator flow carried out so far, apart from the case of 
small flow rates, a simplifying but not strictly justified 
assumption was used. The contribution of a given spin 
generator to the a-effect was considered independent of 
the neighboring spin generators and in that sense deter- 
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FIG. 1: The dynamo module (after Q). The signs + and 
- indicate that the fluid moves in the positive or negative z— 
direction, respectively, in a given spin generator. R — 0.85 m, 
H = 0.71m, a = 0.21m. 




mined under the condition that all its surroundings are 
conducting fluid at rest |2^, p3j, p5[ . An analogous 
assumption was used in calculations of the effect of the 
Lorentz force on the fluid flow rates in the channels of 
the spin generators It remained to be clarified 

which errors result from these assumptions. 

The main purpose of this paper is therefore the calcula- 
tion of the a-coefficient and a related coefficient as well 
as the quantities determining the effect of the Lorentz 
force on the fluid flow rates for an array of spin gen- 
erators, taking into account the so far ignored mutual 
influences of the spin generators. In Section |0] the mod- 
ified Roberts dynamo problem with the spin generator 
flow is formulated. In Section III the numerical method 



used for solving this problem and the related problems 
occurring in the following sections are discussed. Section 
IV presents in particular results concerning the excita- 
tion condition for the dynamo with spin generator flow. 
In Section |^ various aspects of a mean-field theory of 
the dynamo experiment are explained and results for the 
mean electromotive force due to the spin generator flow 
are given. Section VI deals with the effect of the Lorentz 



force on the flow rates in the channels of the spin gener- 
ators. Finally in Section VII some consequences of our 



FIG. 2: The Roberts flow pattern. The flow directions cor- 
respond to the situation in the dynamo module if the co- 
ordinate system coincides with that in Fig. |l|. 




findings for the understanding of the experimental results 
are summarized. 



FIG. 3: The spin generator flow pattern. As for the flow 
directions the remark given with Fig. |^ applies. The fluid 
outside the cylindrical regions where flow directions are indi- 
cated is at rest. There are no walls between the cells. 



Independent of the recent comprehensive accounts of 
the mean-field approach to the Karlsruhe dynamo ex- 
periment |2^, |23|, ^4), this paper may serve as an intro- 
duction to the basic idea of the experiment. However, 
we do not strive to repeat all important issues discussed 
in those papers, but we mainly want to deliver the two 
supplements mentioned above. 



II. FORMULATION OF THE DYNAMO 
PROBLEM 

Let us first formulate the analogue of the Roberts dy- 
namo problem for the spin generator flow. We consider a 
magnetic field B in an infinitely extended homogeneous 
electrically conducting fluid, which is governed by the 
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induction equation, 

?/V 2 B + V x (it x B) -d t B = 0, V-B = 0, (1) 

where r\ is the magnetic diffusivity of the fluid and u 
its velocity. The fluid is considered incompressible, so 
V • u = 0. Referring to the Cartesian co-ordinate system 
(x, y, z) mentioned above we focus our attention on the 
cell < x, y < a and introduce there cylindrical co- 
ordinates (r, if, z) such that the axis r = coincides with 
x = y = a/2. We define then the fluid velocity u in this 
cell by 



everywhere 
it z = — u for < r < 7"i 

u v = —lot, u z = —(h/2ir)u) for r± < r < T2 



u r = 
u v = 0, 



(2) 



u v = 0, 



u z = 



for r > r%, 



where u and oj are constants, r% and r 2 are the radius 
of the central channel and the outer radius of the helical 
channel, respectively, and h is the pitch of the helical 
channel. The coupling between u v and u z in ri < r < r 2 
considers the constraint on the flow resulting from the 
helicoidal walls of the helical channel. The velocity u 
in all space follows from the continuation of velocity in 
the considered cell in the way indicated in Fig. [| i.e. 
with changes of the flow directions from each cell to the 
adjacent ones so that the total pattern is again periodic 
in x and y with the period length 2a and independent of 
z. 

We characterize the magnitudes of the fluid flow 
through the central and helical channels of a spin gen- 
erator by the volumetric flow rates Vc and Vh given by 



Vc = Trrju, Vk = ^(r| -r\)hu). 



(3) 



We may measure them in units of arj, so we introduce 
the dimensionless flow rates Vc and Vh, 



Vc = Vc/ar) , Vh = Va/ar] . 



(4) 



We further define magnetic Reynolds numbers i? m c and 
i? m H for the two channels by i? m c = ur\ji) and R m n = 
ujr2(r2 — T\)lr\. Thus we have Vc = {nri/a)R m c and 

V H = [(ri +r2)h/2ar 2 ]RmH- 

In view of the application of the results for the consid- 
ered dynamo problem to the experimental device we men- 
tion here the numerical values for the radius R and the 
height H of the dynamo module, the lengths a, h, r\ and 
T2 characterizing a spin generator and the magnetic diffu- 
sivity 77 of the fluid: R = 0.85 m, H = 0.71 m, a = 0.21 m, 
h = 0.19 m, n/a = 0.25, r 2 /a = 0.5, T) = 0.1m 2 /s. 
(More precisely, the values of R and H apply to the "ho- 
mogeneous part" of the dynamo module, i.e. the part 
without connections between different spin generators. 
The value of rj is slightly higher than that for sodium 
at 120°C, considering the effective reduction of the mag- 
netic diffusivity by the steel walls of the channels.) The 
given data imply ar] = 75.6 m 3 /s. Furthermore we have 



Vc = 0.785i? m c and Vh = l-357i? m H, so Vc and Vh are 
in fact magnetic Reynolds numbers. Concerning devia- 
tions from the rigid-body motion of the fluid assumed 
here and the role of turbulence we refer to the more com- 
prehensive representations [^2| p3| |. 

We are interested in dynamo action of the fluid motion, 
so we are interested in growing solutions B of ([!]) with 
the velocity u defined by (Q) and the explanations given 
with them. According to some modification of Cowling's 
anti-dynamo theorem growing solutions B independent 
of z are impossible; cf. ^6|. We restrict our attention to 
solutions of the form 



B = Re[B(x, y, t) exp(ifcz)] 



(5) 



where B is a complex periodic vector field which has 
again a period length 2a in x and y, and k a non- 
vanishing real constant. In this case we may consider 
equations (|l|) in the period interval —a < x,y < a only 
and adopt periodic boundary conditions. (Solutions B 
with larger period lengths, as were investigated for the 
Roberts flow [^7], |29| , seem to be well possible but are 
not considered here.) 

If we put B(x,y,t) = B(x,y) exp(pi) with a parame- 
ter p, for which we have to admit complex values, equa- 
tion ([l]) together with the boundary conditions pose an 
eigenvalue problem with p being the eigenvalue parame- 
ter. Clearly p depends on Vc, Vh and k. The condition 
Re(p) = defines for each given k a neutral line, i.e. a 
line of marginal stability, in the VcVh -diagram, which 
separates the region of Vc and Vh in which growing B 
are impossible from that where they are possible. 



III. THE NUMERICAL METHOD 

In view of the numerical solution of the induction equa- 
tion (Fy) we express B by a vector potential A, 



B = V x A 



(6) 



Inserting this in (Q) and choosing V • A properly we may 
conclude that 

v V 2 A + uxB-d t A = 0. (7) 
Analogous to (JH)) we put 

A(x, y, z, t) = Re[A(x, y, t) exp(ifcz)] . (8) 
Then we have 

B = VxA + ifexA, (9) 

where k = ke with e being the unit vector in z-direction, 
and 



?y(V 2 - k 2 )A + u x B - d t A = . 



(10) 



With a solution A we can calculate B according to (||) 
and finally B according to (0). 
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In the sense explained above we consider (|10|) only 
in the period unit —a < x,y < a and adopt peri- 
odic boundary conditions. When replacing A(x,y,t) by 
A(x, y) exp(pi) we arrive again at an eigenvalue problem 
with p as eigenvalue parameter. 

Let us, for example, assume that p is real and consider 
the steady case, p — 0. We may then consider, e.g., 
Vc as eigenvalue parameter while Vh and k are given. 
Modifying the equation resulting from (Q) by an artificial 
quenching of Vc and following up the evolution of A, the 
wanted steady solutions of the original equation ([!(]) and 
thus the relations between Vc and Vh for given k and 
p = can be found. 

For the numerical computations a grid-point scheme 
was used. They were carried out on a two-dimensional 
mesh typically with 60 x 60 or 120 x 120 points, and some 
of the results were checked with 240 x 240 points. The 
x and y-derivatives were calculated using sixth order ex- 
plicitly finite differences, and the equations were stepped 
forward in time using a third order Runge-Kutta scheme. 
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FIG. 4: Neutral lines describing steady dynamo states in the 
Vc Vh- plane for various values of k. 



V. THE MEAN-FIELD APPROACH 



IV. THE EXCITATION CONDITION OF THE 
DYNAMO 



Using the described numerical method solutions B of 
the dynamo problem posed by ([!]) , (g) and (||) have been 
determined. As in the case of the Roberts flow |2^| the 
most easily excitable solutions are non-oscillatory, which 
corresponds to real p, and possess a contribution inde- 
pendent of x and y. 

Fig. H shows the neutral lines in the VcVh -diagram 
for several values of the dimensionless quantity n defined 
by k — ak. In view of the Karlsruhe experiment the 
case deserves special interest in which a "half wave" of 
B fits just to the height H of the dynamo module, so k, — 
ira/H = 0.929. The neutral line for this case can provide 
us a very rough estimate of the excitation condition of 
the Karlsruhe dynamo. However, this estimate neither 
takes into account the finite radial extend of the dynamo 
module nor realistic conditions at its plane boundaries. 
Let us consider, e.g., the values of Vh necessary for self- 
excitation in the experimental device for given Vc- The 
values of Vh obtained in the experiment as well as those 
found by direct numerical simulations are by a factor 
of about 2 higher than the values concluded from the 
neutral curve for k = 0.9; see e.g. Fig. 4 in Refs. 0] and 
§, Fig. 2 in Ref. @ or Fig. 3 in Ref. ||. The tendency 
of the variation of Vh with Vc is however well predicted. 
(The influence of the finite radial extend of the dynamo 
module on t he excitation condition will be discussed in 
Section VD. It makes the mentioned factor of about 2 



plausible.) 



The Karlsruhe dynamo experiment has been widely 
discussed in the framework of the mean-field dynamo 
theory; see e.g. Q. Let us first discuss a few aspects of 
the traditional mean-field approach applied to spatially 
periodic flows and then a slight modification of this ap- 
proach, which possesses in one respect a higher degree 
of generality. We always assume that the magnetic flux 
density B is governed by the induction equation ([j]) and 
the fluid velocity u is specified to be either a Roberts 
flow or the spin generator flow as defined above. 



A. The traditional approach 

For each given field F we define a mean field F by tak- 
ing an average over an area corresponding to the cross- 
section of four cells in the xy-plane, 



F(x,y,z) = ^ j I F(.r -<..;/-:- <<///. (11) 



We note that the applicability of the Reynolds averag- 
ing rules, which we use in the following, requires that F 
varies only weakly over distances a in x or y-direction. 
(The following applies also with a definition of F using 
averages over an area corresponding to two cells only pl| , 
but we do not want to consider this possibility here.) 

We split the magnetic flux density B and the fluid 
velocity u into mean fields B and u and remaining fields 
B' and u', 

B = B + B', u = u + u'. (12) 
Clearly we have u = 0, and therefore u = u 1 . 
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Taking the average of equations (P we see that B has 
to obey 

r)V 2 B + V x £ -d t B = 0, V-B = 0, (13) 
where £, defined by 



£ = u x B' , 



(14) 



is a mean electromotive force due to the fluid motion. 

The determination of £ for a given u requires the 
knowledge of B' . Combining equations (]l|) and ( [ilf ) we 
easily arrive at 

r]V 2 B' + V x (u x B'Y - d t B' = -Vx(mxB), 

V-B' = 0, (15) 



where (uxB')' = uxB' —u x B' . We conclude from this 
that B' is, apart from initial and boundary conditions, 
determined by u and B and is linear in B. We assume 
here that B' vanishes if B does so (and will comment on 
this below) . Thus £ too can be understood as a quantity 
determined by u and B only and being linear and ho- 
mogeneous in B. Of course, £ at a given point in space 
and time depends not simply on u and B in this point 
but also on their behavior in some neighborhood of this 
point. 

We adopt the assumption often used in mean-field dy- 
namo theory that B varies only weakly in space and time 
so that B and its first spatial derivatives in this point 
are sufficient to define the behavior of B in the relevant 
neighborhood. Then £ can be represented in the form 



£i = 



b ijk dBj/dxk , 



(16) 



where the tensors ay and byk are averaged quantities 
determined by u. We use here and in the following the 
notation xi = x, X2 = y, X3 = z and adopt the summa- 
tion convention. Of course, the neglect of contributions 
to £ with higher order spatial derivatives or with time 
deri vativ es of B (which is in one respect relaxed in Sec- 
tion VB) remains to be checked in all applications. 

The specific properties of the considered flow patterns 
allow us to reduce the form of £ given by (|l^) to a more 
specific one. Due to our definition of averages and the 
periodicity of the flow patterns in x and y, and its inde- 
pendence of z, the tensors ay and bijk are independent 
of x, y and z. Clearly a 90° rotation of the flow pattern 
about the z-axis as well as a shift by a length a along 
the x or y-axes change only the sign of u so that simul- 
taneous rotation and shift leave u unchanged. This is 
sufficient to conclude that ay and byk are axisymmetric 
tensors with respect to the z-axis. So ay and bijk con- 
tain no other tensorial construction elements than the 
Kronecker tensor <5; TO , the Levi-Civita tensor e; mn and 
the unit vector e in z-direction. The independence of 
the flow pattern of z requires that ay and bijk are i n ~ 
variant under the change of the sign of e. Finally it can 
be concluded on the basis of (111) that £ has to vanish if 



B is a homogeneous field in z-direction, which leads to 
<Z33 = 0. With the specification of ay and bijk according 
to these requirements relation (16) turns into 



£ = -a±(B-(e-B)e) 

-P±V xB (fli - P ± )(e • (V x B)) e 
-ftex(V(e.5) + (e-V)B), (17) 

where the coefficients aj_, /3m and (3^ are averaged 
quantities determined by u and independent of x, y and 
z. The term with a± describes an a-effect, which is 
extremely anisotropic. It is able to drive electric currents 
in the x and y-directions, but not in the z-direction. The 
terms with /3± and give rise to the introduction of a 
mean-field diffusivity different from the original magnetic 
diffusivity of the fluid and again anisotropic. In contrast 
to them the remaining term with /3 3 is not connected 
with VxB but with the symmetric part of the gradient 
tensor of B and can therefore not be interpreted in the 
sense of a mean-field diffusivity. 

In the case of the Roberts flow the coefficient a± has 
been determined for arbitrary flow rates, and coefficients 
like /3_l, /3|| and /3 3 for small flow rates fig |o|, ||, H> ||- 
As for the spin generator flow only results for a± have 
been given so far jl9|, ^(| |25| . 

For the determination of a±_ it is sufficient to consider 
equation ( |l5| ) for B' with B specified to be homogeneous 

field. In this case, which implies V x (u x B 1 ) = 0, this 
equation turns into 

rjSI 2 B' + (B' ■ V)u - (u ■ V)B' - d t B' = - (B ■ V)u , 



V-J3' = 0. 



(18) 



We may again consider B' like B as independent of z. 
Let us put B = B ± + B u and u = uj_ + uii with 
B ± = B — (e • B')e and B m = (e ■ B')e, and u± and 
defined analogously. Then we find 

r?V 2 .B' ± + {B' ± ■ V)w± - (ul • V)B' ± - d t B' ± 

= -(B-V)u x , 
ilV 2 B\\ - (u x ■ V)B'n - d t B'\\ (19) 

= -((S + B' x )-V)tt,|. 

We further put Uj_ = u±uj_ and uu — u\\Ui\, where u± 
and W|| are factors independent of x and y characterizing 
the magnitudes of u± and «|| , and u± and M|| fields which 
are normalized in some way. Clearly B' ± is independent 
of uii, and B || linear in uii. The x and y-components 

of u x JB', from which a± can be concluded, are sums of 
products of components of ity and B' ± and of u± and 
B ii. Thus a± must depend in a homogeneous and linear 
way on u\\ whereas the dependence on u± is in general 
more complex. This can be observed from the results for 
the Roberts flow. In view of the spin generator flow we 
split My into two parts, ttyi and 1*1121 °f which the first 
one is non-zero in the central channel and the second 
one in the helical channel only. We further introduce 
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the corresponding quantities u\\ i and u\i 2 - We may then 
conclude that a± is linear but no longer homogeneous 
in uii i. Since uu i is proportional to Vc we find that a± 
is linear but not homogeneous in Vc whereas it shows a 
more complex dependence on Vh- 

For small flow rates we may neglect the terms with u 
on the left-hand side of equation ( |l8| ) . This corresponds 
to the second-order correlation approximation often used 
in mean-field dynamo theory. Then the solutions B' and 
further a± can be calculated analytically. Starting from 
the result found in this way for the spin generator flow 
fj"9"| , |2(i| , p5| and using the above findings we conclude that 
the general form of ax reads 



a 2 hr) 



[VcMVa/h^ + ^VnMVu/hr])] (20) 



with two functions 0c and 0h satisfying 0c (0) = 
0h(O) = 1. Note that the argument Vn/hr] is equal to 
(a/h)Vn, which is in turn equal to uj(r 2 +ri)(r 2 —r\)j2r\. 
Consequently it is just some kind of magnetic Reynolds 
number for the rotational motion of the fluid in a he- 
lical channel. The functions 0c and 0h have been 
calculated analytically under a simplifying assumption 
p0| , p5[ , which, however, proved not to be strictly cor- 
rect. We will give rigorous results for a± and f or 0c and 
0h in a more general context below in Section V C . 

As announced we make now a comment on the as- 
sumption that B' vanishes if B does so. Investigations 
with the Roberts dynamo problem have revealed that 
non-decaying solutions B of the induction equation (|l|) 
whose average over a cell vanishes are well possible [ B8j . 
They coincide with non-decaying solutions B' of the 
equation (|l^) in the case B = 0. These solutions are, 
however, always less easily excitable than solutions with 
non-vanishing averages over a cell. They are therefore 
without interest in the discussion of the excitation condi- 
tion for mean magnetic fields B. In that sense the above 
assumption is, although not generally true, at least in 
the case of the Roberts flow acceptable for our purposes. 
Presumably this applies also for the spin generator flow. 

In view of the next Section we assume for a moment 
that B does not depend on x and y but only on z. In that 
case we have VxB = ex [V(e-£?) + (e-V)-B] = exdB/dz 
and therefore ( |l7| ) turns into 



£ = - a± (B- (e-B)e) 



fie x dB/dz, 

fi = fix + fis ■ (21) 



Interestingly enough, here the difference in the characters 
of the fix and /33-terms in (|l7]) is no longer visible. While 
there are reasons to assume that the coefficients fix and 
fiu , which can be interpreted in the sense of a mean-field 
diffusivity, are never negative, this is no longer true for fi^ 
and therefore also not for fi. The results for the Roberts 
flow show indeed explicitly that fi can take also negative 
values ITg|,|22l M. 



B. A modified approach 

We now modify the mean-field approach discussed so 
far in view of the case in which B does not depend on x 
and y but may have an arbitrary dependence on z. All 
quantities like B, B, B or 5, which depend on z, are 
represented as Fourier integrals according to 



F(x,y,z,t)— / F(x,y,k,t)exp(ikz)dk . 



(22) 



The corresponding representation of B clearly includes 
the ansatz (|^). B depends on x, y, k and £, but B 
and £ depend only on k and t. The requirement that 
F(x, y, z, t) is real leads to F*(x, y, k, t) = F(x, y, —k, t). 
Relations of this kind apply to JB, B, B' and £. 

Equations ([ll]) to ( |i"5| ) remain valid whereas ( fl6|), (JTj 
and (21) have to be modified. Clearly (p"3|) and (|14[) are 



equivalent to 

■qk 2 ~B~ -ikx£ + d t B~ = . 

and 



e • B = 0. 



£ = u x B 



(23) 



(24) 



Instead of (15) we have 

7?(V 2 - k 2 )B' + (V + ife) x (u x B)' - d t B 
= -(V + ife) x (u xS), (V + ifc)-i?' = G, (25) 



where (u x B )' = u x B u x B . 

Assuming again that £ is linear and homogeneous in 

B we conclude that the same applies to £ and B, too. 
Therefore we now have 

£ l (k,t) = a lJ {k)~B J (k,t), (26) 

where djj is a complex tensor determined by the fluid 

flow. Analogous to £ and B it has to satisfy a*j(k) — 
Q,j(-i). From the symmetry properties of the it-field 
we conclude again that the connection between £ and 
B remains its form if both are simultaneously subject 
to a 90° rotation about the z-axis, i.e. relation ( p6f ) re- 
mains unchanged under such a rotation of £ and B. This 
means that the tensor is axisymmetric with respect 
to the axis defined by k. The general form of that is 
compatible with a*j(k) = aij(-k) is given by 

dy(fc) = ai(\k\)Sij + a 2 (\k\)kikj + ia 3 (|fc|)e y 7&z (27) 

with real a\, a 2 and 03. Together with ( ^6| ) this leads to 
£z = (0.1 + CL2k 2 )B z . On the other hand £ z is equal to 
the average of u x B' y — u y B' x , and we may conclude from 

( p5| ) that B' x and B' y are independent of B z . This in turn 
implies a\ +a 2 k 2 — 0. We note the final result for &ij(k) 
in the form 



&ij(k) = -&x{k)(Sij - e t e-j) + \fi(k)eijih 



(28) 
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with two real quantities a± and 0, which are even func- 
tions of k. 

From (Eq) and rt28) we obtain 



£(k) = -a±{k)(B~ - (e ■ B)e) - 0(k)k x B . (29) 
Together with (|2|) this leads to 
£{z,t) = 

a± (k) [B(k, t) - (e • B(fc, t))e] exp(ifcz)dfc 
- e x y (3(k)B~(k, t) exp(ifcz)dfc . (30) 
This in turn is equivalent to 

/ ax(0[S(z + C,0-(c-5(« + C,*))e]dC 



with 



-^ex^ / >\!U(:-,./)<k 



ax(C) = / <5x(&) exp(ifc£)dfc . 
f3{C) = { (3(k)exp(ik()dk. 



(31) 



(32) 



Note that both a± and (3 are even in £. 

Let us now expand &ij(k) as given by (28) in a Taylor 
series and truncate it after the second term, 



&ij(k) = -a±(0)(Sij - e^j) + ife/3(0)e ij ie i 



(33) 



The corresponding expansion of £ as given by (30) reads 
£ = -ol(0)(B- (e-B)e) -/3(0)e X dB/dz. (34) 

Comparing this with relation (||l]) of the preceding sec- 
tion we find 



«x = «±(Q), /9 = j8(0). 



(35) 



Returning again to arbitrary k we define for later pur- 
poses a function d(k) by 



d(fc) = Q ± (fc) + k(3{k) . 



(36) 



If a{k) is given, we may determine dx(fc) and /3(fc) ac- 
cording to 



ox(fc) = ^ [a(k) + &{—k)] , 



Moreover, we have 



-[&(k) - a(-k)] . 



ax = fi(0), /?=^(0). 



(37) 



(38) 



C. The parameters defining a— effect etc. 

In view of the determination of the quantities ox(&) 
and /3(fc), which includes that of ox and (3, we note that 
relations like (|l]) or (^) connecting £ with B or £ with 

-B apply, apart from the explicitly mentioned restrictions, 
for arbitrary B. Thus we may take these quantities from 
calculations carried out for specific B. 



Using the method described in Section III we have nu- 
merically determined steady solutions of equation ( |l5| ) 

for B with given Vq , Vh, k and a specific B of Beltrami 

type satisfying e x dB/dz = kB. With these solutions 

we have then calculated the quantity £ ■ B , which, ac- 
cording to (|j)|) , has to satisfy 



SB = 



\B\ 



(39) 



with a defined by (j36|). From the values of a and their 
dependence on k obtained in this way a±(k), /3(fc), a± 
and (3 have been determined. 

In mean-field models of the Karlsruhe device in the 
sense of the traditional approach explained in Section 
VA the coefficient a±_ occurs in the dimcnsionless com- 



bination C — a±R/r), with R being the radius of the 
dynamo module, and the influence of (3 can be discussed 
in terms of (3 = /3/rj. We generalize the definitions of C 
and (3 by putting 



(40) 



Now C and f3 show a dependence on fc, which we express 
by one ontc = ak. 

Thinking first of the traditional approach we consider 
C with k — 0. Figure H shows contours of C in the Vc^h - 
diagram, Fig. |^ the functions <fic and </>h, from which a± 
and thus C can be calculated. These results deviate for 
large Vk significantly from those determined with the 
simplifying assumption mentioned above, according to 
which the mutual influence of the spin generators was 
ignored (2(], P5| . In the region of Vc and Vh which is 
of interest for the experiment, say < Vc, Vn < 2, the 
values of C, 4>c and 0h for given Vc and Vh are somewhat 
larger than those obtained with that assumption. One 
reason for that might be that in the case of an array 
of spin generators, compared to a single one in a fluid at 
rest, the rotational motion in a helical channel expels less 
magnetic flux into regions without fluid motion, where it 
can not contribute to the a-effect. Remarkably, in the 
region < Vc , Vh < 2 our result for C agrees very well 
with one derived under the assumption of a Roberts flow 

HE!- 

Figure ^ exhibits contours of (3 for n = in the VcVr- 
diagram. We already pointed out that (3 can take neg- 
ative values. Here we see that (3 becomes negative for 
sufficiently large values of Vc and Vh- Although this 
happens somewhat beyond the region of interest for the 
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FIG. 5: Contours of C for k = 0. 



experiment it suggests that inside this region the posi- 
tive values of (3 may be small. The diffusion term in the 
mean-field induction equation is proportional to rj(l+/3). 
In the investigated region of Vc and Vh this quantity 
proved always to be positive. 

Let us now proceed to C and (3 for k ^ 0. As already 
mentioned, in view of the experimental device it seems 
reasonable to put k — ira/H = 0.929. Analogous to 
Fig. [| which applies to k = 0, Fig. || shows contours of 
C for k = 0.9. We see that C for given Vc and Vh is 
slightly higher in the latter case. The results for (3 are 
virtually indistinguishable for both cases. 



D. The excitation condition in mean— field models 
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FIG. 6: The functions <j>c and <j>n calculated for an array of 
spin generators. For comparison the results of the approxi- 
mation considering single spin generators (i.e. ignoring their 
mutual influences) are also given. 



We consider first again the traditio nal a pproach to 
mean-field theory explained in Section V A . Equation 
( |l3| ) for B together with relation (|l~7|) for £ allows the 
solutions 



B = £?o( cos(fcz), =F sin(fcz), 0} 
p = -(-n + (3)k 2 ±a±k, 



exp(pt) , 



(41) 



where Bq is an arbitrary constant. We refer here again 
to Cartesian co-ordinates and consider k as a positive 
parameter. For these solutions we have V x B = ±kB, 
i.e. they are of Beltrami type. This implies that there are 
no mean electric currents in the z-direction. The solu- 
tion that corresponds to the upper signs can grow if a±_ 
is sufficiently large. The condition of marginal stability 
reads a±_ = (77 + j3)k or, what is the same, 



C= {l+/3)kR, 



(42) 



where C and /3 have to be interpreted as the values for 
k = 0. If we relate this to the dynamo module and put 
k = tt/H we have 



C= (l+0)vR/H. 



(43) 



V c [m 3 /h] 
100 200 300 




Vr 



FIG. 7: Contours of /3 for k = (calculated as the limit 
K -> 0). 
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FIG. 8: Contours of C for k = 0.9. 

Note that the factor R in the conditions (^) and (Q) 
results from the definition of C only. In fact they are 
independent of R. 

Proceeding to the modified approach to the mean-field 
theory and replacing relation (17^ for £ by ( |30| ) we find 
formally the same result. However, a± and (3 have to be 
replaced by a± and /3, and C and (3 in ( f42|) a nd (Q) have 
to be taken for k = ak. The conditional 2j) interpreted 
in this sense defines neutral lines in the VcVh -diagram 
which have to agree exactly with those shown in Fig. ^. 
Likewise, the condition ( |43"| ) defines the special neutral 
line with n = Tta/ H. 

One of the shortcomings of estimates of the self- 
excitation condition of the experimental device based on 
the solutions of the induction equation used in Section 



IV or, equivalently, on a relation like (|43|), consists in 
ignoring the finite radial extent of the dynamo module. 
We point out another solution of equation (|l3|) for B 7 
which has been used for an estimate of the self-excitation 
condition of the experimental device considering its fi- 
nite radial extent [£0[ plf . For the sake of simplic- 
ity we assume that £ is given by equation ( p~7| ) with 
(3± = j3« =^3=0. We refer to a new cylindrical co- 
ordinate system (r, tp, z) adjusted to the dynamo module 
so that r = coincides with its axis and z = with its 
midplane. The solution we have in mind reads 

B = B — , — — (r&) exp(pt , 

oz a_L r or I 



9 = Jo(qr) cos(fcz) , 



(44) 



where q and fc are constants and Jo is the zero-order 
Bessel function of the first kind. This solution is axi- 
symmetric with respect to the z-axis. It has further the 
property that the normal components of V x B vanish 
both on the cylindrical surfaces qr = z u , where z v denotes 
the zeros of Jo, and on the planes kz = (I + l/2)n with 



integer I. We identify the region inside the smallest of 
these cylindrical surfaces between two neighboring planes 
of that kind with the dynamo module, so we put q = 
Zi/R, where z\ is the smallest positive zero of Jo, and k = 
Tt/H. Then there are no electric currents penetrating the 
surface of the dynamo module. The condition of marginal 
stability for the so specified solution reads 



C = n(R/H)[l + ( Zl H/nR) 2 ] . 



(45) 



In the limit H/R — > this agrees with (|43|) if we put 
(3 = 0. For finite H/R, however, C is now always larger 
than the value given by ([ilf ) with (3 = 0. This can easily 
be understood considering that there is now an additional 
dissipation of the magnetic field due to its radial gradient. 
C as function of H/R has a minimum at H/R = ~k / z\. 
The dynamo module was designed so that H/R has just 
this value. In this case we have 



C = 2nR/H . 



(46) 



In other words, the real radial extent of the dynamo 
module enlarges the requirements for C, compared to 
the case of infinite extent, by a factor 2. As can be 
seen from Fig. |^, in the region of Vc an d Va m which 
experimental investigations have been carried out, say 
1.3 < Vc, Vk < 1-6, this enlargement of C means that if, 
e.g., Vc is given, Vh grows by a factor between 2.5 and 
3.5. We recall here the deviation of the experimental re- 
sults from the estimate of the self-excitation condition 
given in Section IV on the basis of Fig. |J, which just cor- 
responds to (E3j) . In the light of these explanations con- 



cerning the influence of the radial extent of the dynamo 
module this deviation is quite plausible. It is actually 
rather small, which indicates that our reasoning despite 
a number of neglected effects does not underestimate the 
requirements for self-excitation. 

We note that also the result ([l6|) is not a completely 
satisfying estimate of the self-excitation condition of the 
experimental device. Apart from the fact that it does not 
consider realistic boundary conditions for the dynamo 
module it is based on an axisymmetric solution of the 
equation for B. Several investigations have however re- 
vealed that a non-axisymmetric solution is slightly easier 
to excite than axisymmetric ones 19, [2l| ^2| ^3). The 
influence of the (3± and /?3-terms of £ can no longer be 
expressed by (3, and there is also an influence of the /3y- 
term. All these influences increase the marginal values of 

c m. 



VI. THE EFFECT OF THE LORENTZ FORCE 
ON THE FLOW RATES 

In the theory of the experiment, equations determining 
the fluid flow rates in the loops containing the central 
channels and in those containing helical channels have 
been derived from the balance of the kinetic energy in 
these loops. The rate of change of the kinetic energy in 
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a loop is given by the work done by the pumps against 
the hydraulic resistance and the Lorentz force. For the 
work done by the Lorentz force averaged over a central 
or a helical channel we write (it • f)V, where (• • •) means 
the average over this channel, V its volume and / the 
Lorentz force per unit volume, 



/ = ^(VxB)xB, 



(47) 



with /i being the magnetic permeability of free space. 

We use again B = B + B' . For all results reported 
here we have assumed that B is a homogeneous field 
and, correspondingly, B' is also independent of z so that 
equations (1^) apply. Then also / is independent of z 
and (• • •) may simply be interpreted as an average over 
the section of the channel with the xy-plane. 

We have calculated the quantities (u- f)c and (it • /}h 
for a central and a helical channel analytically in two 
different approximations |2^, |2^| . In approximation (i) 
all contributions to / of higher than first order in Vc or 
Vh were neglected so that it applies to small Vc and Vh 
only. In approximation (ii) arbitrary Vc and Vh were 
admitted, but as in an earlier calculation of the a effect 
only a single spin generator surrounded by conducting 
medium at rest was considered, i.e. any influence of the 
neighboring spin generators was ignored. We represent 
the results of both approximations in the form 



(u-f)c 
(u-/) H 



a (V c 



27c \sc 



5iV>c(Vc,Vk) : 



— (—) siVBiw.r,,). mm 



Here a is the electric conductivity of the fluid, 7c and 
7h are given by 



7c 



1 



7h 



(n + r 2 ) 2 + (V^) 2 

2(r 2 +r 2 ) + (W 



(49) 



sq and sh are the cross-sections of the central and helical 
channels, and B± is the magnetic flux density perpendic- 
ular to the axis of the spin generator, i.e. to the z-axis. In 
approximation (i) we have ipe — ipn = 1. In approxima- 
tion (ii) ipe and i/>h are functions of Vc and Vh satisfying 
^ c (Vc,0) = 1 for Vc ^ and ipn{V C ,0) = 1 for all Vc, 
varying only slightly with Vq and decaying with growing 
Vh; see also Fig. The factors ipe and ipn in the rela- 
tions ( ff8| ) for (it • /)c and {u ■ /)h describe the reduction 
of the Lorentz force by the magnetic flux expulsion out 
of the moving fluid by its azimuthal motion. 

We may conclude from the relevant equations that 
( u ' /)c and {u ■ /}h can again be represented in the 
form ( |48| ) if the complete array of spin generators and 
arbitrary Vc and Vh are taken into account. Only the 
dependences of ipe and V'h on Vc and Vh changes. 

Before giving detailed results we make a general state- 
ment on these dependences. As in the considerations in 
the paragraph containing (|l8| ) we may again introduce 
the quantities B\, JSV uj_, uu and use (pi). With the 



same reasoning as applied there we find that for the spin 
generator flow B ± is independent of Vc and B 11 linear 
in Vc- We further express f ± and defined analogous 
to B'± and £?'||, according to (|47j ) by the components of 
B' j_ and B » , their derivatives and the components of B. 
In this way we find that {u ■ f)c is a sum of two terms, 
one proportional to Vc and the other proportional to Vq. 



Consequently, ipe has the form Tp^\Vn) + ip c i; (Vh)V c 1 

with V c 0) (°) = !■ We further find that (u ■ f) h is a 
sum of three terms, one independent of Vc and the oth- 
ers proportional to Vc and Vq, and ipn has the form 

V-h 0) (Vh) + iI$(Vh)V c + ^h 2) (Vh)-V c 2 with ^ 0) (0) = 1. 
This can be seen explicitly from the calculations in the 
approximation (ii) mentioned above, in which, by the 

/ ( 2 ) n 

way, Vh = °- 

We have calculated tpc and i/'h numerically on the basis 
of equations (|l5| ) using the method described in Section 
III. The result is shown in Fig. ^. Instead of the com- 
plete array of spin generators we have also considered 
an array in which fluid motion occurs only in one out 
of 4 x 4 spin generators. The numerical result obtained 
for this case agrees very well with the analytical result of 
approximation (ii) shown in Fig. |^. 

For a complete array of spin generators the factors tpc 
and ipn in the relations (Eq) are generally larger com- 
pared to approximation (ii). In other words, the Lorentz 
force is less strongly reduced by the azimuthal motion of 
the fluid. This can be understood by considering that 
less magnetic flux can be pushed into regions without 
fluid motion. 



VII. CONCLUSIONS 

We have first dealt with a modified Roberts dynamo 
problem with a flow pattern resembling that in the Karls- 
ruhe dynamo module. Based on numerical solutions of 
this problem a self-excitation condition was found. Since 
in these calculations neither the finite radial extent of 
the dynamo module nor realistic boundary conditions at 
its plane boundaries were taken into account this self- 
excitation condition deviates markedly from that for the 
experimental device. 

A mean-field approach to the modified Roberts dy- 
namo problem is presented. Two slightly different treat- 
ments are considered, assuming as usual only weak vari- 
ations of the mean magnetic field in space, or admitting 
arbitrary variations in the ^-direction. The coefficient 
a_i_ describing the a-effect and a coefficient (3 connected 
with derivatives of the mean magnetic field are calcu- 
lated for arbitrary fluid flow rates. The result for aj_ 
corrects earlier results obtained in an approximation that 
ignores the mutual influences of the spin generators [p5| . 
It leads to a much better agreement of the calculated 
self-excitation condition with the experimental results 
p2| , |2l| . We note in passing that in the case of small 
flow rates our result, although calculated for rigid-body 
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FIG. 9: The dependence of tpu and tpc on Vh for Vc = 1 
and Vc = 2 for an array of spin generators. For comparison 
the results of approximation (ii) , which considers a single spin 
generator only, are also given. 



motions only, applies also for more general flow profiles 



p2| , . The result for j3 suggests that the enlargement of 
the effective magnetic diffusivity by the fluid motion can 
be partially compensated by another effect of this mo- 
tion. The same has been observed in investigations with 
the Roberts flow [^9). This could be one of the reasons 
why the results calculated under idealizing assumptions, 
in particular ignoring the effect of the mean field diffu- 
sivity, deviate only little from the experimental results 




In the framework of the mean-field approach we have 
also given an estimate of the excitation condition which 
considers the finite radial extent of the dynamo module. 
It shows that the real extent enhances the critical value 
of C, which is a dimcnsionlcss measure of by a factor 
2. In other words, if in the region of V*c and Vh, in which 
experimental investigations have been carried out, Vc is 
fixed, Vh has to be larger by a factor between 2.5 and 3.5. 
If the excitation condition is corrected in this way it does 
not underestimate the requirements for self-excitation. 

We have also calculated the effect of the Lorentz force 
on the fluid flow rates in the channels of a spin genera- 
tor. Again our result corrects a former one obtained in 
the approximation already mentioned which ignores the 
mutual influences of the spin generators J22|, ^3|. The 
braking effect of the Lorentz force proves to be stronger 
than predicted by the former calculations. This means in 
particular that estimates of the saturation field strengths 
given so far [^2[ |24| have to be corrected by factors be- 
tween 0.8 and 0.9; for more details see Note added in 
proof in (24). 
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